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Abstract. The origin of astrophysical magnetic fields observed in galaxies and clusters of galaxies 
is still unclear. One possibility is that primordial magnetic fields generated in the early Universe 
provide seeds that grow through compression and turbulence during structure formation. A cosmo- 
logical magnetic field present prior to recombination would produce substantial matter clustering at 
intermediate/small scales, on top of the standard inflationary power spectrum. In this work we study 
the effect of this alteration on one particular cosmological observable, cosmic shear. We adopt the 
semi-analytic halo model in order to describe the non-linear clustering of matter, and feed it with the 
altered mass variance induced by primordial magnetic fields. We find that the convergence power 
spectrum is, as expected, substantially enhanced at intermediate/small angular scales, with the exact 
amplitude of the enhancement depending on the magnitude and power-law index of the magnetic 
field power spectrum. Specifically, for a fixed amplitude, the effect of magnetic fields is larger for 
larger spectral indices. We use the predicted statistical errors for a future wide-field cosmic shear 
survey, on the model of the ESA Cosmic Vision mission Euclid, in order to forecast constraints on 
the amplitude of primordial magnetic fields as a function of the spectral index. We find that the am- 
plitude will be constrained at the level of ~ 0.1 nG for ub ~ —3, and at the level of ~ 10~ 7 nG 
for us ~ 3. The latter is at the same level of lower bounds coming from the secondary emission 
of gamma-ray sources, implying that for high spectral indices Euclid will certainly be able to detect 
primordial magnetic fields, if they exist. The present study shows how large-scale structure surveys 
can be used for both understanding the origins of astrophysical magnetic fields and shedding new 
light on the physics of the pre-recombination Universe. 
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1 Introduction 

Magnetic fields are ubiquitous in the Universe. They are found, with varying intensities, in all kinds of 
astrophysical objects, ranging from planets (B ~ 1 G) to massive clusters of galaxies (B ~ 1 fiG, see 
[1] for a review). The origin of astrophysical magnetic fields is still debated, and two main scenarios 
have been devised. According to the first, magnetic fields form inside stars through mechanism 
such as the Biermann battery [2], they are then amplified due to dynamo mechanisms [3-5], and 
are finally diffused at larger scales thanks to supernova winds and Active Galactic Nuclei (AGN) 
feedback. The second scenario postulates instead that magnetic fields of very low intensity were 
generated in the early Universe, during inflation [6, 7] or phase transitions [8] (for recent reviews 
see [9, 10]), and were subsequently amplified during gas infall within dark matter halos due to the 
conservation of the magnetic flux and the onset of turbulence [11]. One interesting aspect of this latter 
scenario is that magnetic fields present in the early Universe are capable of generating and amplifying 
density perturbations in the baryonic fluid through the Lorentz force. Thanks to the gravitational 
coupling, these density perturbations would propagate to the dark matter fluid as well. In fact, for 
some time it was believed that such magnetic fields could be responsible for the generation of the 
seed density fluctuations for structure formation [12], until it was realized that the resulting power 
spectrum would be totally different from the observed one [13, 14]. Still, if Primordial Magnetic 
Fields (PMFs henceforth) existed, they would have altered the distribution of density fluctuations 
laid down during inflation. 

The presence of magnetic fields in the early Universe has been considered in the literature as 
early as in [15-17]. Subsequently, the impacts of such fields on the primordial nucleosynthesis and 
the formation of the Cosmic Microwave Background (CMB henceforth) have been thoroughly studied 
[18], and constraints on their strength given by the observed elemental abundance [19-23] and CMB 
temperature power spectrum [24-30] have been put (see also [31] for forecasted PMF bounds from 
Planck). These constraints are in general compatible with a magnetic field strength, extrapolated at 
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present time, that is < InG. More recent works studied the impact of PMFs of this level on several 
aspects of structure formation, such as the formation of voids [32], the redshift-space correlation 
function of galaxies [33, 34], cosmic reionization [35, 36], and 21-cm fluctuations [37]. 

The only direct way of measuring the strength of magnetic fields on cosmological scales is the 
Faraday rotation that these fields induce on the polarization direction of high-redshift sources, such 
as galaxies and quasars. An early example of such a measurement can be found in Vallee [38], who 
reported an upper limit on the strength of cosmological magnetic fields of ~ 0.3 nG (0.023/Slb,o^ 2 )> 
assuming that all baryons are completely ionized in the low-redshift Universe (note that the author 
actually quotes a lower value, which is however based on an electron number density substantially 
higher than today's commonly accepted value). Earlier works (e.g. [39, 40]) as well as the subsequent 
review by Kronberg [41] seemed however to point toward values somewhat higher than this. In fact, 
[42] found an upper limit on the magnetic mean-field amplitude of ~ 0.5 nG (0.023/f2b 5 o^ 2 ) on 
scales ~ 10 — 50 hr x Mpc, while [43] reported findings varying from ~ 1 — 10 nG, depending on 
the coherence length they considered. More recently [44] found an upper limit of ~ 0.8 fiG for the 
intergalactic magnetic field strength on a scale of ~ 1 /i -1 Mpc. 

Several more exotic possibilities for probing cosmological magnetic fields exist. For instance, 
in [45] the authors have proposed to use the CMB itself as a source of polarized radiation in order 
to measure the Faraday rotation due to PMFs. At around the same time, Plaga [46] put forward the 
possibility to use the secondary emission of gamma-ray sources such as gamma-ray bursts (see also 
[47]) or blazars in order to put constraints on the strength of intergalactic magnetic fields. As it turns 
out, the PMF amplitude should be at least of ~ 10 -8 — 10 -6 nG in order to explain this kind of 
observations [48-50]. 

In the present paper we consider a novel way to investigate PMFs. As mentioned above, cos- 
mological magnetic fields present before recombination are expected to alter the large-scale matter 
distribution. Specifically, PMFs introduce substantial power on top of the inflationary one at interme- 
diate/small scales. Cosmic shear, that is the weak gravitational deflection of light by the Large-Scale 
Structure (LSS) should be ideally suited to probe those scales. Here we compute the cosmic shear 
power spectrum produced in cosmological models with PMFs having a variety of spectral slopes. 
Then, using the specifications of future wide-field weak lensing surveys on the model of Euclid 1 
[51], a recently accepted ESA Cosmic Vision program mission, we estimate the statistical error on 
the power spectrum and hence derive forecasted constraints on the amplitude of PMFs. A similar 
study has been recently performed in [52], which we discuss in Section 5. 

The rest of the paper is organized as follows. In Section 2 we review the impact of PMFs 
on the linear stage of structure formation, with particular care to the growth factor of magnetically- 
induced matter density fluctuations and their linear power spectrum. In Section 3 we show how the 
modifications at the linear level can be mapped to the non-linear level, describing the computation 
of the mass function, halo bias, and eventually the non-linear matter power spectrum in presence of 
PMFs. In Section 4 we turn attention to the convergence power spectrum and how it gets modified 
by PMFs, thus giving forecasted constraints on the amplitude of PMFs as a function of their spectral 
index. In Section 5 we summarize our conclusions. 

In this work we adopted the cosmological parameters suggested by the latest WMAP data re- 
lease in conjunction with Baryon Acoustic Oscillation (BAO) and Hubble constant measurements 
[53]. This implies a matter density parameter of f2 m> o = 0.272, a dark energy density parame- 
ter of Qa,o = 1 — ^m,o> a baryon density parameter of $7b,o = 0.046, and a Hubble constant of 
Ho = h 100 km s _1 Mpc -1 , with h = 0.704. The amplitude of primordial curvature perturbations 

'http://www.euclid-ec.org/ 
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is set to = 2.44 x 10 -9 , corresponding to eg = 0.809 in absence of PMFs, while the spectral 
index is given by n = 0.963. Note that PMFs have impacts on the CMB power spectrum as well. 
These display themselves both on scales smaller than those probed by WMAP, due to vortical per- 
turbations [27, 54], and on very large angular scales, due to anisotropic magnetic stresses [55, 56]. 
While the former do not affect the WMAP bounds on cosmological parameters [10], the latter, being 
at the percent level, might alter the inferred amplitude and slope of the inflationary curvature power 
spectrum. Awaiting for more detailed studies on this topic, in this work we chose to stick with the 
WMAP cosmological parameters. If not specified otherwise, we shall assume natural units such that 
the speed of light is c = 1 throughout. 



2 Linear structure growth 
2.1 Magnetic fields 

PMFs affect the growth of density perturbations in the baryon fluid. Since density fluctuations in 
the dark matter component are sourced also by those in the baryonic component, magnetic fields 
indirectly affect the growth of dark matter perturbations as well. Let us indicate with B(x,t) the 
magnetic field that is immersed in the baryon fluid. For the sake of completeness, it is instructive 
to review the equations describing the evolution of small comoving perturbations for a magnetized 
fluid in an expanding cosmological background [12]. While the continuity and Poisson equations are 
identical to the standard ones, the Euler equation displays a significant difference, 

—L±l + H (t)v(x,t) = r±-±J- + S(x,t) , (2.1) 

at a(t) 

where 

q/ ,x [W xB(x,t)]xB(x,t) 

In the previous equations we have neglected the role of pressure gradients, since these are irrele- 
vant on scales much larger than the thermal Jeans length. It appears that the peculiar velocity of 
baryon density perturbations is not only sourced by potential perturbations, but also by magnetic 
fields through the Lorentz force. 

Magnetic fields in turn obey the standard solenoidal condition V • B(x, t) = 0, and also they 
are fueled by the velocity fields via magnetic induction, through a backreaction mechanism that for a 
fluid with infinite conductivity reads 

d r 9, s , m V x \v(x,t) x a 2 (t)B(x,t)] 

- {a*(t)B(*, t)] = 1 > a(t) H ■ ("J 

On scales much larger than the magnetic Jeans length (see below) the backreaction given by the right- 
hand-side of the previous equation can be neglected at first order, thus leading to the constancy of 
a 2 (t)B(x, t), that is a 2 (t)B(x, t) = B(x, to) = Bq(x). From this, we note the interesting fact that 
the magnetic energy density, 

u ,_ {B 2 (x,t)) (B 2 (x)) p Bfl 

PB{t) = ~^r~ = = 5*w ( } 

(where angular brackets represent ensemble averages), evolves in the same way as the radiation en- 
ergy density in an expanding Universe [13]. 
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One important ingredient in order to estimate the impact of PMFs on the growth of cosmic 
structures is their correlation function. For a statistically homogeneous and isotropic magnetic field 
the correlation function in Fourier space can be written as [13, 57] 

'Biik^Bjik'^)) = ^-S D (k - k') Uj - ^) P B (k,t) . (2.5) 
It is often assumed that the power spectrum P B (k, t) of PMFs is a power law, 

P B (k,t) = A B (t)k nB , (2.6) 

in a given range of scales, k G [/c m in 5 ^max]. and vanishes otherwise. For the sake of completeness, in 
this work we considered a wide range of spectral indices. However we bear in mind that constraints 
from the gravitational wave background have been claimed to prefer n B < [58], although these 
constraints are conditional on independent determinations of the effective number of neutrino species. 
Also, joint constraints coming from the CMB and the linear matter distribution favor n B < [29, 59], 
however such constraints strongly depend on the priors assumed in order to break the degeneracy 
between the magnetic amplitude and spectral index. 

The ultraviolet cutoff scale k max is often interpreted as the scale at which magnetic fields are 
dissipated due to radiative viscosity around recombination. Also, the primordial magnetic fields are 
usually assumed to be correlated up to very large scales, so that, effectively, k m - m = 0. Before going 
with more detail in the specific value of /c max adopted in this work, we note that, analogously to 
what happens for density fluctuations, the normalization of the magnetic field power spectrum can 
be written in terms of the variance of the magnetic fields smoothed on a certain comoving scale A. 
Assuming a sharp Fourier-space smoothing function of width k = 2ir/\ < k m3iX we obtain 



Note that some authors (e.g., [56]) do not divide by the factor 2 in Eq. (2.5), so that for a fixed 
spectral amplitude their magnetic field variance is twice as large. This is a point to keep in mind 
when comparing constraints on the amplitude of PMFs obtained in different works. The result in 
the previous equation is valid only if n B > —3, which is a common assumption in PMF studies. It 
follows that a B (X,t) is the magnetic comoving mean-field amplitude across the scale A at cosmic 
time t. 

Now we want to give some more details about the ultraviolet cutoff fc max - According to [14, 60], 
the scale at which magnetic fields suddenly dissipate can be estimated as (see also [56]) 

-^max _ a B (^max)t r ) f tx dt 



fc ma X ( 27r ) 2 4vrp 7ir CT T Jo a 2 {t)n e (t) ' 
In the previous equation p J>T is the radiation energy density at recombination, <tt is the Thomson 
scattering cross section, and n e (t) is the number density of free electrons before the recombination 
time i r . We recall that for a fully ionized plasma with standard primordial abundances n e (t) = 
0.52 ph{t) I fim p , with p = 0.59. Considering a generic length scale A < A max we get 

ALx _4(A,*r) / A \ 3+nB pm p 



) "^^f 1 ^- (2 " 9) 
J 0.52 p bi0 Jo 



(27r) 2 47rp 7jr cr T \Xi 

It is safe to assume that before recombination the Universe behaves as an Einstein-de Sitter universe, 
so that the integral in the previous equation can be easily written as 
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Jo 



a(t)dt= / -™-= J* (2.10) 

Jo H \ a ) 5iioV S 'm,0 

Also, the radiation energy density and the magnetic field variance have the same temporal evolution, 
so that the ultraviolet cutoff scale can be estimated as 

& = -^ o(A) A3 " " mp 4 /2 ■ (2.1D 

5 p^ a T if o yn^0.52p bi0 r 

By placing the relevant numerical values where necessary, in particular a T = 1/(1 + z x ) — 9.083 x 
10 -4 , we obtain 



A, 



5+n s 



A 



~ 1.319 x 10 -3 h 2 



Q"b,o(A) 
InG 



A 



flb,oh 



-1/2 



l/r^Mpc ) V 0.023 J V 0.135 

(2.12) 

As an example, if we assume a typical value for the present-day magnetic field dispersion o\b,o(A) = 
0.1 nG for A = 1 /i _1 Mpc and h = 0.704 we obtain A max ~ 18.7 /i _1 kpc if n B = -2, A max ~ 
91.8 fo -1 kpc if n B = 0, and A max ~ 182 /i _1 kpc if n B = 2. These are the values of the ultraviolet 
cutoff that shall be used in the remainder of this paper. We note that A max is a rather small scale, hence 
any cosmologically relevant scale over which we wish to smooth the magnetic field itself would be 
larger than this. 

Another important length scale for magnetized fluids is the magnetic Jeans length, A#- This is 
the comoving transition scale below which perturbations in the magnetized fluid stop growing, due 
to the fact that magnetic pressure gradients counteract the gravitational pull [14]. This length scale 
has been derived for the first time in [13] for a fully baryonic Universe. Later, the generalization to 
cosmological models dominated by dark matter was provided by [35]. The magnetic Jeans length 
reads implicitly as 



a(t)X B (t) 



2 a B (X B (t),t) 2 a B (X,t) 



A 



X B (t) 



(3+n B )/2 



(2.13) 



5 V G Pm(t)pb(t) 5 ^/Gp m (t)p h (t) 

If the dispersion of the magnetic field evolves as a B (X,t) oc a~ 2 {t), this implies that the magnetic 
Jeans length is independent of time, so that we can define X B (t) = X B: o = X B . It follows that 



~ 0.2298 h 2 

It is easily seen that the magnetic Jeans scale is much larger than the ultraviolet cutoff scale. 

There is one further potentially important length scale for the problem at hand, that is the stan- 
dard thermal Jeans length for baryons. As mentioned above, the magnetically induced density fluctu- 
ations in dark matter are sourced by their baryonic counterparts. On scales smaller than the thermal 
Jeans length the latter are washed out by thermal pressure, hence the former cannot grow either. Nor- 
mally the thermal Jeans scale is much smaller than the magnetic one, and in fact it has been ignored 
in the vast majority of works on PMFs. However, it has recently been argued (see [35, 61, 62]) that 
the small-scale dissipation of PMFs due to ambipolar diffusion and decaying turbulence provides an 
extra source of energy that is injected in the baryonic fluid. This additional energy would increase the 
temperature of the inter-galactic medium by several orders of magnitude at low redshifts, and hence 



VB,o{X) 

1 nG 



J. 



A 



l/i-iMpc 



0.023 



0.135 



(2.14) 
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increase substantially the thermal Jeans scale. The computations of [61] show that the thermal Jeans 
scale can indeed become larger than the magnetic one if one considers the strengthening of magnetic 
fields due to baryonic matter compression in collapsing structures. At the linear level however, the 
increase in temperature would be much lower than that, and hence it would be safe to assume the 
thermal Jeans scale to be always smaller than or comparable to the magnetic one. This is the assump- 
tion that we adopted in the rest of the paper. Plus, the extra heating due to matter (and hence magnetic 
field lines) compression occurring inside the first non-linear structures is arguably compensated by 
the cooling due to the much more abundant H2 [61]. 

2.2 Time evolution of density fluctuations 

Since dark matter is not directly coupled to magnetic fields, the evolution of linear dark matter density 
perturbations is given by the standard equation 

5.(x,t) + 2H(t)S,(x,t) - 4irG[p.(t)5.(x,t) + p h (t)5 h (x,t)} = 0. (2.15) 

In other words, density fluctuations in both baryons, S\,(x,t), and dark matter, 5,(x,t), source the 
growth of the latter due to gravitational instability. On the other hand, baryons are coupled with 
magnetic fields, hence the evolution of their density fluctuations follow a different equation [13], 

5 h (x,t) + 2H(t)S h (x,t) - 4TrG[p.(t)5.(x,t) + Ph (t)5 h (x,t)} = V ' S ^ . {2 A6) 

a[t) 

In this case, the magnetic fields also provide an additional force that contributes to the linear growth 
of density perturbations. By following, e.g., [33] we define the density fluctuation of the total matter 
fluid as 

Pm(t)6 m (x, t) = p.(t)5.(x, t) + p h (t)5 h {x, t) , (2.17) 

where p m = p,+p-^. It easily follows that 5 m (x, t) = f.5,(x, t) + f^,5^(x, t), where /. = fi.,o/fim,o 
and / b = O bj oMn,o- B Y replacing S.(x,t) = 5 m (x,t)/f. - 5 h {x,t)f h /f. into Eq. (2.15) and 
making use of Eq. (2.16) we obtain 

5 m (x,t) + 2H(t)5 m (x,t) - A7rG Pm (t)5 m (x,t) = f h V ' . (2.18) 

a(t) 

Eq. (2.18) is what one needs to solve in order to determine the time evolution of density fluc- 
tuations in the presence of magnetic fields, and hence in order to find out the growth factor. This 
equation can be further simplified by recalling that on the scales A 3> Xb of our interest the mag- 
netic field evolves simply as B(x, t) = Bq{x) / a 2 (t), so that the forcing term of Eq. (2.18) can be 
rewritten as /^V • So(x)/a 3 (t). Now, the homogeneous part of the previous equation is the standard 
equation for the evolution of density fluctuations on large scales. Its solution can be written as a 
linear combination of the usual growing and decaying modes of structure growth, D + (t) and -D-(i). 
In order to obtain the most general solution of Eq. (2.18) we need to find a particular solution of the 
same equation and sum it to the general solution of the homogeneous equation. It follows that the 
solution is 

5 m (x, t) = T+(x)D+(t) + F_(x)D_(t) + / b V • S (x)M(t) ~ T + (x)D + (t) + / b V • S (x)M(t) . 

(2.19) 
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Figure 1. The growth factors relevant for the present work, normalized to unity at z = 0. The black solid line 
represents the growth factor of standard inflationary density perturbations, while the red dashed line shows the 
growth factor of density perturbations induced by PMFs. 

The function M(t) represents the temporal evolution of the matter density fluctuations induced by 
PMFs, and it is a solution of the differential equation 

M(t) + 2H(t)M(t) - AnGp m (t)M(t) = ~4rr ■ (2.20) 

Such a solution has to be found through numerical integration for a ACDM Universe. At very early 
times however the Universe is well represented by an Einstein-de Sitter model, for which the previous 
equation admits an analytic solution for t > t r (before recombination magnetically induced density 
perturbations cannot grow because of the tight coupling between baryons and photons), that is [13] 

M(t) = t 2 

We used this solution in order to set up the initial conditions for the numerical integration of the 
differential equation. It is interesting to note that the growing part of the analytic solution above is 
identical to the standard growth factor in an Einstein-de Sitter universe. By numerically integrating 
the differential equation we verified that this the case also for the more general ACDM case. 

In Figure 1 we compare the standard growth factor D + (z) for inflationary density perturbations 
with the growth factor of density perturbations induced by PMFs, M(z). At low redshift the two 
growth factors coincide, as is the case for the Einstein-de Sitter universe. In fact, according to Eq. 
(2.21), M(z) oc (1 + z)^ 1 once the non-growing part has been dissipated. This is due to the fact that, 
once the density fluctuations have been generated by magnetic fields, they grow due to gravitational 
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Figure 2. The power spectrum of density fluctuations induced by PMFs. Each line style and color refer to a 
different magnetic spectral index, as labeled, while the spectral amplitude is set by <tb,o(1 /i~ 1 Mpc) = 0.1 nG. 
The black solid line shows the power spectrum of primordial density fluctuations for reference. 

instability as the primordial ones. The main difference between the two function stays at high redshift, 
and in particular in the fact that, while magnetically induced density perturbations cannot grow before 
recombination (Log(l + z r ) ~ 3), the standard density perturbations have been growing for a while 
by z r , so that the decaying mode has already dissipated. As a consequence, we have that M (z) — > 
as z — > z T . Note that in the redshift range we shall be interested in this work (z < 5), the two growth 
factors can be effectively considered as identical. 

2.3 Matter power spectrum 

In order to evaluate the linear matter power spectrum in presence of PMFs we make the assumption 
that primordial density fluctuations are uncorrelated from the density fluctuations induced by mag- 
netic fields. This is an assumption that is often employed in the literature, although the presence or 
absence of such correlations depends on the mechanism responsible for the generation of PMFs, over 
which a consensus has not been reached yet (see [63] for a discussion). If however that is the case, 
then it is easily seen that the total matter power spectrum is 

P m (k, t) = D%(t)P(k) + M\t)U(k) , (2.22) 

where P(k) is the standard matter power spectrum, while TT(fc) is the power spectrum of matter den- 
sity fluctuations induced by magnetic fields. In writing down the previous equation we assumed that 
the decaying mode has dissipated, so that the growth of standard density perturbations is effectively 
given only by the growing mode. It has been shown in the literature [13] that the power spectrum of 
magnetically induced fluctuations can be written as 
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n(fc) = /b „ / °° dg / d p B(q)PB[a(k,q,ri] ^ + ^ _ 2 + ^3 5 3, 
(47rpb,o) -/o V-i a 2 (k,q,n) 

(2.23) 

where a(k,q,fi) = ^Jk 2 + q 2 — 1kq\i and /1 is the cosine of the angle between k and q, that is 
k ■ q = kqfi. We note that, since the magnetic field power spectrum is affected by the ultraviolet 
cutoff at /c max , the two integrals in the equation above actually do not cover the entire reported 
integration range. We are now in a position to estimate the total linear matter power spectrum and its 
evolution with cosmic time. 

In Figure 2 we show the contribution of magnetically induced density fluctuations to the total 
lineal - matter power spectrum, for different values of the magnetic spectral index ranging from ng = 
— 2 to riB = 2. Here and in what follows, unless specified otherwise, we normalized the magnetic 
field power spectrum such that o\b,o(1 fa Mpc) = 0.1 nG. The generic effect of PMFs is to sharply 
increase the linear matter power spectrum at scales k ~ 1 — 10 fa Mpc -1 . On scales k > fc max 
the impact of PMFs abruptly vanishes, due to the absence of magnetic field correlations above that 
wavenumber. In practice the magnetically induced linear power spectrum should be cut off at the 
magnetic Jeans scale ks, which is significantly lower than /c max (see Eqs. 2.12 and 2.14). In Figure 
2 we did not include this cut-off, for demonstration reasons. However for the computation of the 
mass variance (Section 3 below) and all non-linear quantities that follow we did introduce the Jeans 
length cut-off. We finally note the fact that at large scales the slope of the power spectrum induced 
by PMFs is shallower for tib = —2 than for larger spectral indices. This is expected, as in [33] the 
authors shown that for k/k mSLX <C 1 the matter power spectrum is n(/c) oc k 4 if ng > —3/2, while 
U(k) oc k 2nB+r if n B < -3/2. In the specific case of ns = -2, we get U(k) oc k 3 . 



3 Non-linear structures 



3.1 Variance 

An essential ingredient in order to estimate the mass function and bias of dark matter halos is the 
mass variance. In cosmological models with PMFs the variance has to be computed from the total 
lineal - matter power spectrum given in Eq. (2.22), that is 

1 r+00 

a 2 {R,z) = ^J k 2 dkP m (k : z)W 2 {k,R) , (3.1) 

where W(k, R) is the Fourier transform of a smoothing function, that we assumed to be a top-hat 
of radius R in real space. In Figure 3 we show the rms (the square root of the variance) of the 
total density fluctuation field as a function of the mass corresponding to a given smoothing scale, at 
z = 0. Due to the small-scale bump in the matter power spectrum produced by PMFs, the variance 
approaches the inflationary value for high masses, while it tends to grow much faster than that and 
reach much larger values when moving to smaller and smaller masses. 

For the calculation of the variance, and for all applications that follow, we truncated the magneti- 
cally induced linear power spectrum at the magnetic Jeans wavenumber ks- This is an approximation, 
since a complete treatment taking into account the backreaction of velocity fields for k ~ ks [13] 
results in a growth factor that depends on scale and indeed vanishes for k > ks, but not with a sharp 
transition. Yet, this is a good approximation and has been used regularly in the literature [35, 49]. 
Because of this cut-off, in models with nonvanishing PMFs the variance (and hence the mis) flattens 
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Figure 3. The mean amplitude of the total matter density field smoothed on a scale corresponding to mass 
m, at z = 0. Different line styles refer to different slopes for the magnetic field power spectrum, while the 
amplitude is set by 0\b,o(1 /i _1 Mpc) = 0.1 nG. The black solid line shows the result obtained by ignoring 
PMFs. 

out for masses smaller than that corresponding to the magnetic Jeans length. This magnetic Jeans 
mass can be estimated as 



' M B 
M(A) 

where 



~ 1.214 x 10" 2 h 6 



Q"b,o(A) 
1 nG 



A 



1 fe" 1 Mpc 



-6 



0.023 



0.135 



(3.2) 



M(A) = 



^m,0i? 2 A 3 „ 3161 x 10 H fc -l M . 



2G 







A 



n 



m.O 



0.272 



(3.3) 



1 fc^Mpc, 

is the average mass contained within the radius A. For ctb,o(1 fa _1 Mpc) = 0.1 nG and the slope 
values considered here the magnetic Jeans mass ranges between Mb ~ 10 9 — 10 10 h~ 1 M & . 

By looking at the curves that show the effect of PMFs in Figure 3, we note some oscillations in 
the mean density amplitude at low mass values. The reason for these are the oscillations present in 
the Fourier transform of the smoothing function, combined with the sharp cutoff of the magnetically 
induced matter power spectrum at the magnetic Jeans scale. As a matter of fact we verified that if the 
top-hat window function is replaced by a Gaussian (as many authors in the literature do, e.g. [36]), 
the oscillations disappear. However, it should be kept in mind that virtually all the mass function 
theory is based on a top-hat window function, hence we chose to stick with this choice. 
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Figure 4. The cumulative mass function for dark matter halos at z = 0, as a function of mass. The black 
solid line refers to the standard case with no PMFs, while the other lines refer to magnetic fields with different 
spectral slopes, as labeled. In all cases the magnetic field power spectrum is normalized by requiring that 
o-B,o(l h^Mpc) =0.1nG. 

A few consideration are in order about the behavior of the mass variance described above. 
Since for high masses the variance is basically unaffected by the presence of PMFs, we expect the 
abundance of massive galaxy clusters not to be substantially modified, and if so, only in cases in 
which the amplitude or spectral index of the magnetic field power spectrum are large enough. On the 
other hand, one might argue that the abundance of low-mass dark matter halos (m < 10 12 /i -1 M Q ) 
should be largely increased. However, there is another effect that needs to be taken onto account. 
The mass function is exponentially sensitive to the mass variance only for the rarest objects at any 
given redshift. This means that, since at z = low-mass objects are not in the exponential tail of the 
mass function, the impact of PMFs on their abundance is also expected to be relatively modest, and 
dictated only by the larger slope of the rms. 

3.2 Halo Mass Function and Bias 

We now turn attention to the estimates of dark matter halo mass function and bias. A proper treat- 
ment of these issues would need fully numerical cosmological simulations run with the modified 
initial conditions induced by PMFs. We plan to explore this topic in a future investigation, since this 
load of work is beyond the scope for the present paper. As a viable alternative, we adopted prescrip- 
tions based on semi-analytic considerations. Specifically, we employed the Sheth & Tormen mass 
function [64] and the Sheth, Mo, & Tormen linear bias [65] prescriptions. Since these formalisms are 
based at least partially on physical motivations (specifically, the ellipsoidal collapse model), we can 
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Figure 5. The linear bias factor for dark matter halos at 
to the standard case with no PMFs, while the other lines 
labeled. In all cases the magnetic field power spectrum 
nG. 



z = 0, as a function of mass. The black solid line refers 
refer to magnetic fields with different spectral slopes, as 
is normalized by requiring that <tb, o(l h~ 1 Mpc) = 0.1 



argue them to remain acceptably valid in cosmologies including PMFs, as long as the mass variance 
a 2 (m, z) is computed properly as described in Section 3.1. 

In Figure 4 we show the mass function obtained with the Sheth & Tormen prescription at z = 
for PMFs whose amplitude is set by (Jb,o(1 h~ 1 Mpc) = 0.1 nG and for different spectral slopes. We 
plot the cumulative mass function, defined as 

N(m,z) = / d£n(£,z) . (3.4) 

J m 

As can be seen, the mass function is basically not modified at very high masses. This is due to the 
fact that at those masses the variance is also unaffected, because PMFs do not alter the linear matter 
power spectrum at large scales. The halo number counts get however increased by a factor of ~ 5 
at intermediate masses, i.e., at m ~ 10 13 — 10 14 /i _1 M Q for rig = 2 and lower masses for lower 
magnetic spectral indices. At very low masses the effect of PMFs is again reduced (and eventually 
even reversed), due to the fact that those masses are not at the tip of the mass function at z = 0. 
Therefore their abundance is given by the slope of the variance as a function of mass (which flattens 
out below the magnetic Jeans mass) rather than the absolute value of the variance itself. 

In Figure 5 we report the linear bias of dark matter halos, computed according to the Sheth, Mo, 
& Tormen prescription. In a fashion reminiscent of the mass function, the bias is almost unchanged 
at low and high masses, while it is substantially reduced (by up to 60%) at m ~ 10 13 — 10 14 h^ 1 M Q 
for ?ib = 2. Also, the lower the magnetic spectral index, the lower the mass where the maximum 
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reduction is observed, again similarly to the behavior of the mass function. This can be easily inter- 
preted as a consequence of the increase in the variance shown in Figure 3 being shifted at lower and 
lower masses as the magnetic spectral index decreases. Additionally, the fact that the halo bias drops 
in the same regime where the mass function increases is self-consistent, as more abundant objects 
have a lower clustering strength. 

3.3 Non-linear Power Spectrum 

In order to study cosmic shear, a model for the non-linear clustering of matter is needed. In this work 
we employed the halo model [66-68], a physically motivated framework that allows to write down 
the fully non-linear matter power spectrum as the sum of two terms, Pnl (k, z) = P\ (k, z) +P2 (k, z), 
where 



+00 r */-__ i.\t2 



Pi(k,z)= / dmn(m,z) 







p(m, z, k) 

Pm,0 



and 



P 2 (k,z) = P m (k,z) 



' , , u p(m,z,ky 
am n{m, z) o{m, z)~ 



(3.5) 



(3.6) 



/O Pm,0 

The first term defines the contribution from dark matter particle pairs belonging to the same halo, 
while the second term includes the clustering of different halos. Let us briefly review the ingredients 
of the model, while for details on the implementation we refer to [69]. The function b(m, z) repre- 
sents the linear halo bias, while n(m, z) is the halo mass function. We already described how both 
were modeled in Section 3.2. The function p(m, z, k) is the Fourier transform of the average dark 
matter density profile, p(m, z, r), while Pm{k, z) is the linear matter power spectrum, discussed in 
Section 2.3. 

As the presence of PMFs introduces more matter power at small scales, it is likely that the 
internal structures of dark matter halos, and hence their average density profile, are modified with 
respect to the standard cosmological model without magnetic fields. Still, a full account of this mod- 
ification could be obtained only with fully numerical cosmological simulations. Since these are not 
available yet, we assumed throughout that the average dark matter density profile is well represented 
by a NFW function [70]. Moreover, we connected the concentration of the halo to its virial mass by 
adopting the model described in [71]. This is admittedly an uncertainty of our approach, however in 
the subsequent cosmic shear analysis we limited ourselves to large and intermediate scales, so that 
the impact of the very inner structure of dark matter halos is not going to be relevant. 

In Figure 6 we compare the linear matter power spectrum in absence of PMFs with its non- 
linear counterpart obtained via the halo model. We also show the two contributions to the non-linear 
power spectrum, the 1— halo term, dominating at small scales, and the 2— halo term, dominating at 
large scales. One point worth noting is that, as well known, the 1— halo term at large scales tends to 
a constant, namely 

r+00 m 2 

\imPi(k,z)= / dmn(m,z)—^ — . (3.7) 
JO P 2 m ,0 

Hence, the higher the mass function at high masses (note the m 2 weighting), the higher this constant. 
This is going to be relevant for the subsequent discussion. 

In Figure 7 we show how PMFs affect the non-linear matter power spectrum. In this plot the 
amplitude of PMFs is fixed, while results for different spectral indices are shown. As can be seen, the 
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Figure 6. The linear power spectrum and its non-linear counterpart obtained via the halo model at z = 0, as 
labeled. PMFs are ignored here. The two black solid lines represent the 1— halo (dominating at large k) and 
2— halo (dominating at small k) contributions to the non-linear spectrum. 

non-linear clustering of matter increases substantially at intermediate and small scales. One part of 
this increment is due to the corresponding increment that is seen in the linear matter power spectrum, 
shown in Figure 2. Please note that here we applied the sharp cutoff at the magnetic Jeans scale, so 
that the impact of PMFs disappears at smaller wavenumbers compared to Figure 2. Another part is 
more indirect, and shows up at scales larger than those affected by the former. This second part is 
due to the increase in the mass function at intermediate masses displayed in Figure 4, which results 
in an increase of the small-fc plateau of the 1— halo term, as shown in Eq. (3.7). In other words, the 
1— halo term becomes important at larger scales as compared to the case with no PMFs, and this is 
more and more pronounced as the spectral index or the spectral amplitude become larger. 

The physical interpretation of this last fact is that structure formation starts earlier in cosmolo- 
gies with PMFs than in the standard model, due to the enhanced mass variance that allows density 
perturbations to surpass the threshold for collapse at an earlier time. This means that density fluc- 
tuations with a given physical scale enter earlier in the non-linear regime or, turning the argument 
around, fluctuations that are still evolving linearly at a given redshift in the standard cosmology have 
already entered the non-linear evolution phase when PMFs are introduced. This finding is new and 
very important. Up to now the impact of PMFs on galaxy clustering has been ignored because it was 
argued that such an impact would manifest itself only at scales too small to be important. This is true 
for the linear power spectrum, but Figure 7 shows that including non-linear modifications to the mat- 
ter clustering can invalidate this assumption. More specifically, PMFs with large enough amplitude 
and/or spectral index may be able to affect the BAO region. This means that, in principle, a study 
of the BAO could be used to put constraints on cosmological magnetic fields. However this kind of 
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Figure 7. The linear (black dotted line) and non-linear (black solid line) matter power spectra computed in 
absence of PMFs. The colored dashed lines represent the non-linear power spectra computed with PMFs 
having various spectral slopes, as labeled. In the latter cases the PMF amplitude is always set by requiring that 

o-b,o(1 /i -1 Mpc) =0.1nG 

study would require a more detailed modeling of the BAO scales, with either numerical simulations 
or perturbation theory. This is an interesting idea for future investigation. 

4 Cosmic shear 

4.1 Convergence power spectrum 

In order to derive the power spectrum of effective convergence [72] from the fully non-linear matter 
power spectrum we adopted the Limber's approximation, which is valid for scales that are sufficiently 
small for the sky to be considered flat. Jeong & Komatsu [73] have shown that this approximation is 
good at the percent level for angular scales 9 = 2-k jl < 36 degrees. Accordingly, the convergence 
power spectrum reads 



In the previous equation \ = x( z ) i s the comoving distance out to redshift z, the scale factor a{x) 
is normalized such that a(0) = 1, and Jk{x) 1S the comoving angular diameter distance, which in 
general depends on the curvature of the Universe. In this work we considered only flat cosmological 
models, for which /k{x) = X> however we keep the general notation for completeness. 

The function W(x) * s a geometric weight, taking into account the fact that if the deflector is too 
close to the observer or to the sources the lensing strength is suppressed. It reads 
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Figure 8. The convergence power spectrum for the fiducial cosmology without magnetic fields (black solid 
line). The colored dashed lines refer to PMFs with different spectral slopes, as labeled. In all latter cases the 
magnetic field amplitude is set by requiring erg, o(l h~ 1 Mpc) = 0.1 nG. 



where g(z) is the source redshift distribution. All these integrals formally extend up to the comoving 
horizon distance Xh, however in practice the source redshift distribution drops to zero well before 
that distance is reached. In this work we adopted the source redshift distribution prescribed for 
Euclid [51], adopting the configuration that covers 15, 000 square degrees. This is 



where F(a) is the complete Euler gamma function, z* = 0.6374, a = 2, and /3 = 1.5. The corre- 
sponding median source redshift is z m = 0.9. 

In Figure 8 we show the convergence power spectrum computed for the various PMF spectral 
slopes considered in this work. As can be seen, the sharp cutoff corresponding to the magnetic Jeans 
scale is now smoothed out. This is due to the fact that a given angular scale in the sky corresponds to 
different physical scales at different redshifts, so that a given multipole gets contributions from phys- 
ical scales both above and below the cutoff. We remind the reader that the comoving magnetic Jeans 
scale is independent of redshift [13]. Consistently with the effect seen for the matter power spectrum, 
PMFs cause a substantial increase in convergence power at intermediate/small angular scales. For a 
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fixed PMF amplitude, this effect is more marked for larger spectral slopes. This obviously implies 
that constraints on the spectral amplitude will be stronger for larger values of n B . 



4.2 Constraints on PMFs 

In order to obtain bounds on the amplitude of PMFs it is necessary to quantify the statistical error 
on the convergence power spectrum estimated through cosmic shear. If ignoring tomography and the 
non-Gaussian part of weak lensing covariance, that however dominates only at small angular scales, 
such an error can be estimated as [74-77] 



AC K {£) 



A£(2£ + A£)f sky 



C K {£) + ^ 
n 



(4.4) 



In the previous equation / s k y represents the fraction of the sky covered by the weak lensing survey 
at hand, <t 7 is the intrinsic shape contribution to the average galaxy ellipticity, and n is the average 
number density of background sources. The multipole interval A£ is the width of multipole band 
where the power spectrum estimator is averaged upon. We adopted A£ = 1, and for the other 
parameters we chose the Euclid specifications, that are f sky = 0.364, <r 7 = 0.3, and n = 30/arcmin 2 . 

Given this, we can estimate a Ax 2 function by assuming that the measured (and hence true) 
convergence power spectrum is the one obtained in the standard model without PMFs, while the 
model power spectrum is the one obtained by including PMFs with a given amplitude and spectral 
index. Therefore we can write 



A X 2 (A B)0 ,n B ) = 



C K (A Bfi ,n B ;t)-Cj?\£) 
AC^(£) 



(4.5) 



where Ab,o is actually determined by <tb 5 o(1 /i _1 Mpc) and the superscript (0) refers to quantities 
evaluated within the fiducial framework having vanishing PMF amplitude. The choice of the multi- 
pole range where the sum is to be extended is important. The Euclid specifications require £\ = h and 
£2 = 5, 000, however we adopted two more conservative extrema. In particular, we set £\ = 50, in or- 
der to make sure that the Limber's approximation is excellently satisfied, and we adopted £2 = 3, 000 
in order to make sure that the impact of baryonic physics (e.g., [78]), the non-Gaussian part of the 
weak lensing covariance, and the uncertainties on the inner structure of dark matter halos in the pres- 
ence of PMFs can all be safely neglected. Needless to say, if all these uncertainties will be carefully 
modeled in the future, the analysis will be extended to substantially higher multipoles, with resulting 
tighter constraints on the amplitude of PMFs. 

In Figure 9 we show the resulting constraints on the amplitude of PMFs as a function of the 
spectral index, for different confidence levels. We computed these bounds for a few selected values 
of n B and then interpolated over the full range of spectral indices shown. This procedure introduces 
a slight degree of uncertainty, however, given the very fluctuating nature of many of the present 
bounds on the amplitude of PMFs, and because the calculation of the convergence power spectrum 
in these models is quite time-consuming, we believe a full Monte-Carlo exploration of the parameter 
space is not warranted. As one might expect, the constraints are tighter for larger than for smaller 
spectral indices, because the impact of PMFs on the linear matter power spectrum is also larger in 
those cases (see Figure 2). Constraints on the amplitude of PMFs range from ~ 0.1 nG for spectral 
indices close to n B = —3, down to ~ 10™ 7 nG for n B ~ 3. It is worth recalling that observations 
of secondary emission from gamma-ray sources place lower limits on the amplitude of PMFs at the 
level of ~ 10~ 8 — 10~ 6 nG, thus implying that Euclid will necessarily detect PMFs, if they exist 
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Figure 9. Constraints (upper limits) on the PMF amplitude as a function of the spectral index ns, obtained 
from cosmic shear. The different lines refer to different Confidence Levels (CLs). Specifically, the short dashed 
line shows the 99.7% CL, the long dashed line refers to the 95.4% CL, and the solid line displays the 68.3% 
CL. 

and if they have a large spectral index, ng ~ 2 — 3. The saddle point at ng ~ —1.5 is interesting 
and likely due to the change in the slope of the magnetically induced matter power spectrum for that 
specific value of the spectral index. 

5 Summary and conclusions 

In this work we investigated the impact of PMFs on cosmological weak lensing, and assessed the 
power of future wide-field cosmic shear surveys on the model of Euclid for constraining the amplitude 
of these seed magnetic fields. Euclid is a space mission selected for the ESA Cosmic Vision program, 
and it is currently scheduled for launch in 2019. Our main results can be summarized as follows. 

• Owing to the fact that PMFs induce a substantial increase in power at small scales, structure 
formation begins at earlier times if magnetic fields are taken into account. As a consequence, 
fluctuations on scales that at a given redshift are still in the linear regime for the standard 
cosmology with no magnetic fields, are already in the nonlinear evolutionary stage when PMFs 
are introduced. 

• Because of the anticipated non-linear evolution of density fluctuations and the increase in power 
at small scales, PMFs cause a substantial increment in the convergence power spectrum at 
intermediate and small angular scales. The increase is more marked for large spectral indices, 
because the power of matter fluctuations induced by PMFs is larger in those cases. 
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• The former effect allows one to use cosmological weak lensing in order to put constraints 
on the amplitude of PMFs for a given spectral index. In agreement with the previous point, 
constraints are stronger for larger indices: adopting the Euclid cosmic shear specifications 
results in bounds varying from ~ 0.1 nG for rig ~ — 3 down to ~ 1CP 7 nG for ?ib ~ 3. 

Our cosmic shear analysis relies on a series of simplifying assumptions. For instance, only 
the statistical errors on convergence power spectrum estimates were included, while weak lensing 
systematics such as photometric redshift errors and shape measurement uncertainties were neglected. 
Also, the degrading effect of uncertainties on the other cosmological parameters was not taken into 
account. For instance massive neutrinos, acting as a hot dark matter component, can wash away part 
of the small-scale power introduced by PMFs, and hence can somewhat worsen the constraints shown 
in Figure 9. Furthermore, tomography, which has the potential of ameliorating the constraining power 
of weak lensing, was not implemented. These points can be improved upon in future investigations on 
PMFs, while for the sake of the present paper we were more interested in illustrating the approximate 
performance of cosmic shear statistics. 

While this paper was being finalized we became aware of a similar study, performed in [52]. 
There the authors used current cosmic shear data in order to constrain PMFs, by simply summing 
the PMF-induced matter power spectrum to the ordinary non-linear [79] spectrum. The present work 
extends that by [52] by including the impact of PMFs on the non-linear growth of cosmic structures. 
By comparing our results with theirs we can conclude that future cosmic shear surveys will improve 
constraints on PMF amplitude by about one order of magnitude compared to current data. To the 
best of our knowledge ours are the first forecasted constraints that are not based on the CMB, the 
primordial nucleosynthesis, or the linear matter power spectrum. For magnetic spectral indices ng > 
2 the constraints from Euclid will hit the lower limit set by the observations of secondary emission 
by gamma-ray sources, and hence in this case PMFs will be certainly detected, if they exist. It should 
however be recalled that bounds from the CMB and the gravitational wave background seem to favor 
a negative magnetic spectral index, although many uncertainties still remain in this respect. If this is 
the case, then PMFs will be detectable by Euclid only if their amplitude is > 1 pG. 

It is worth recalling that the PMF amplitude and spectral index give important indications about 
the mechanism responsible for their formation, and of the epoch in the early Universe during which 
such formation occurred. For instance, PMF formation occurring during phase transitions is expected 
to produce a bluer spectrum as compared to inflationary formation [10]. Thus, the study of PMFs by 
means of the LSS provides an exciting new window for studying the physics of the early Universe. 
Additionally, it will help to answer the fundamental question of the origin of astrophysical magnetic 
fields observed in galaxies and clusters of galaxies. 
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